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Abstract 

Recent theoretical and experimental studies of the exotic atoms 
with Z = 1 are reviewed. An interplay between the atomic internal 
and external degrees of freedom is essential for a good description of 
the atomic cascade. The perspective of ab initio cascade calculations 
is outlined. 



1 Introduction 

When heavy negative particles ir~ , p, etc.) stop in matter, they usually 
form exotic atoms in highly excited states with principal quantum number 



n ~ ym/m e where m is the reduced mass of the exotic atom and m e is 
the electron mass. The exotic-atom formation is followed by an atomic cas- 
cade consisting of multistep transitions to lower atomic states. For hadronic 
atoms, the atomic cascade is a complete life history because the hadrons 
get absorbed by the nuclei with high probability before reaching the ground 
state. Muonic atoms (where the absorption is weak) de-excite to the ground 
state and engage in various reactions (muon catalyzed fusion, muon trans- 
fer, molecular formation) with initial conditions determined by the atomic 
cascade. In both cases, the atomic cascade can reveal important information 
about the properties of exotic atoms and reactions with atoms in excited 
states. This paper, supplementing the earlier reviews focuses on the 

recent progress in theoretical studies of the atomic cascade in light muonic 
and pionic atoms. 
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Table 1: Cascade processes in exotic atoms with Z — 1 and their energy 
dependence. 



2 Cascade Mechanisms 

A brief summary of the cascade processes in the exotic atoms with Z = 1 
is given in Table |l| The radiative de-excitation and the nuclear absorption 
(in hadronic atoms) do not depend on experimental conditions directly. All 
other processes occur in collisions with surrounding atoms and their rates 
are proportional to the hydrogen density and usually depend on energy. 

At least tree cascade mechanisms are essential for the basic understanding 
of the atomic cascade jjj]: the radiative transitions, the external Auger effect, 
and the Stark mixing. In this paper, the cascade models, which include these 
three mechanisms only, will be called the minimal cascade modelQ (MCM). 

Figure[I| demonstrates the n/-dependence of the total radiative and Auger 
de-excitation rates for muonic hydrogen. The main features of these de- 
excitation mechanisms were discussed in 

The Auger rates calculated in the Born approximation (Fig. |TJh>) are en- 
ergy independent. The eikonal approximation || predicts a rather weak 
energy dependence, with the results being very close to the ones in the Born 
approximation for n < 6 and for a kinetic energy of the order of 1 eV. The 
initial and final state interactions in the Auger transitions were discussed 
in |17||, however, no detailed calculations have been done. 

The Stark mixing corresponds to transitions among the nZ-sublevels with 
the same n. It is a very fast collisional process because the exotic atoms 



1 In the literature, it has also been called the standard cascade model (SCM). 
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Figure 1: The rates of (a) radiative and (b) Auger de-excitation (LHD) in 
muonic hydrogen. 



with Z — 1 are small and electroneutral and have no electrons, so that 
they can easily pass through the regions of the strong electric field inside 
ordinary atoms. When the Stark mixing rate is much larger than all other 
transition rates, the statistical population of the nZ-sublevels is determined 
by the principle of detailed balance. 

The relative importance of the various cascade processes in muonic and 
pionic hydrogen is demonstrated in Figs. |2|a,b. The cascade models of the 
exotic atoms with Z — 1 are listed in Table 0; they are of two types. In 
one group, there are various implementations of the MCM [!],[/], where 
the kinetic energy is assumed to be constant through the whole cascade. 



The other group ||. p9| |31[| consists of detailed kinetics models which take 
the energy evolution during cascade into account by explicit treatment of 
acceleration and deceleration mechanisms. 



3 X-Ray Yields 

Since the rates of the radiative transitions are well known, the competition 
between the radiative and collisional processes can be used for testing the 
collisional de-excitation rates by measuring the X-ray yields at various den- 
sities. The most suitable system for this study is the muonic hydrogen where 
the X-ray yields are not suppressed by absorption during the cascade. An- 
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Figure 2: The effective (statistical average) rates for (a) muonic hydrogen M 
and (b) pionic hydrogen |29j in liquid hydrogen. The effective absorption 
rates for the ir~p defined as in [|IJ and the deceleration rates correspond to 
the kinetic energy E = 2 eV (solid line) and E = 30 eV (dotted line). The 
Coulomb de-excitation rates for E = 2 eV (solid line) and E = 0.04 eV 
(dashed line) are from [20,21 . 
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Table 2: The cascade processes included in various models of the lightest 
exotic atoms in addition to the MCM containing the radiative and Auger 
de-excitation and the Stark mixing. 
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Figure 3: (a) The density dependence of the K a , Kp, and yields in 
muonic hydrogen. The theoretical curves are from ||, the experimental data 
from p2|-|39||. (b) The density dependence of the high-energy components 
(E > 1, 8, 50 eV) in the \ip ground state after the atomic cascade calculated 
in the model M with the Coulomb rates scaled by a factor kc = 0.5 |29 . 



other convenient factor is that the rates of the Auger de-excitation, which is 
the main collisional process, have a weak energy dependence, therefore they 
are not strongly affected by uncertainties in the kinetic energy distribution. 
The main features of the density dependence of the yields of the i^-lines were 
already fairly well explained by the MCM JF|,|8|]. Figure |3]a shows the exper- 
imental data in comparison with the recent calculations [[| that include, in 
addition to the Auger de-excitation, the Coulomb transitions. 

The relative role of the collisional processes can be illustrated by the 
following simplified calculation of the density dependence of the K a /K^ ratio 



using the method suggested in pH| . The yields Y Ka and Y K are given by the 
balance equations: 



A3u a + 0(Agr + x^t mb ) 

xi 



*K a ~ P2 ~ P3 



Yk, = J*%, Af = \U 2 + XU + ^(XtT + &° mb ) (2) 
where p n are the populations of the atomic states n — 2, 3, XJ^ = Ag P _^ ls /3 
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and A^ 2 = + 1)^3/^2p/9 are the effective rates of the radiative transi- 
tions, is the hydrogen density in units of liquid hydrogen density N = 4.3 • 
10 22 cm- 3 (LHD), XfZT and Xg^™* are the rates of the Auger and Coulomb 
transitions normalized to LHD. Equations (0,0) are valid at <p > 0.05 and 
give the following dependence of the ratio Y^/Y^g on the density: 

V \T A( \ Au 9 er I \Coulomb\ / \Coulomb\ 

w, = + aL"" 2 = 079 + u ^ i 1 + IgH (3) 

If one neglects the Coulomb de-excitation, then Eq. (|3|) gives (y# a /Y K _)*' 1 = 
15 at cf) = 1 which is slightly below the corresponding experimental ra- 
tio (Y Ka /Y K0 ) exp = 19.9 ± 2.5 At = 0.08 the difference is larger: 
(Y K JY Kl3 y h = 2.0 and {Y Ka /Y Ka ) exp = 6.4 ± 1.3. With the Coulomb transi- 
tions taken into account the theoretical ratio gets closer to the experimental 
values as discussed in [pl] , |38| . This observation can be considered as evidence 



that some mechanisms in addition to the Auger effect, like the Coulomb tran- 
sitions, are needed for a better description of the X-ray yields. 



4 Kinetic Energy Distribution in Excited States 

The evolution of the kinetic energy distribution during the atomic cascade 
is very important because many collisional processes are energy dependent 
(Table 1). In muonic atoms, it determines the initial energy distribution 



in the ground state, which influences the diffusion of //-atoms pT|-4l3]] and 



muon catalyzed fusion [H4l|4H|. A large fraction of the atoms is not thermal- 
izedQ during the atomic cascade, as it follows from the cascade calculation 
demonstrated in Fig.|3|b. 

The kinetic energy distribution in the atomic cascade can be studied with 
different methods. Direct probes which are model independent are based on 
the measurements of the Doppler broadening of the X-ray lines in fip and 
n~p and the n-TOF spectra in the reaction ir~p — > tt° + n. Given the 
kinetic energy distribution w(E) at the instant of the radiative transition, 
the Doppler broadening of the X-ray line g(Auj) has the form: 



w(E) 



E 



dE 




(4) 



2 The atoms with the kinetic energy much larger than the temperature are called ep- 
ithcrmal and the ones with E> 1 eV — 'hot' or 'highly-energetic'. 
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Figure 4: (a) The kinetic energy distribution w(E) at the instant of the 
3P — > IS transition in muonic hydrogen at 15 bar, (b) the corresponding 
Doppler broadening of the Kp line and (c) the cumulative energy distribution 
W(E). 



where uq is the X-ray energy in the atomic rest frame and M is the mass 
of the atom. A straightforward inversion of the transformation (^) gives the 
cumulative energy distribution W(E): 

W{E) = f / w{E')dE' = 2 / g(x)dx - 2£v / £#(£v / £) . (5) 

JO J 

Figure |3] shows an example of the calculated distributions for the ^-transition 
in fip atoms at 15 bar which features characteristic peaks in the high-energy 
component resulting from the Coulomb de-excitation discussed in Sec. |6|. 
The Doppler broadening of the neutron TOF spectra from the reaction 
(7T~p) — > 7T° + n is related to the kinetic energy distribution in a similar 
way (see |17S-fi8||). A good knowledge of the kinetic energy distribution is 



essential for precise measurements of the {'R~p)is nuclear width [49]. For ex- 
ample, the kinetic energy T = 0.5 eV corresponds to a Doppler broadening of 
the Kp line 5T = 0.1 eV which is about 10% of the nuclear width Vis = 1 eV. 

Indirect methods of probing the kinetic energy distribution rely on models 
of the kinetics. In particular, the first evidence for epithermal muonic atoms 



was found in muon catalyzed fusion (see [1~5] and references therein). Im- 
portant results on the initial kinetic energy distribution in the ground state 
were obtained from the diffusion of fip and fid atoms at low density [^l| -|43 



they allow one to determine the energy distribution in excited states using 
cascade models as discussed in 0,0. Indirect methods exploiting energy 
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dependent processes, like the muon transfer in excited states PB| , |5U| , |5T|1 , were 
used for the estimates of the average energy in excited states, however, they 
are not sensitive to the details of the energy distributions. 



5 Elastic Collisions 

The elastic collisions, where the principal quantum number of the exotic atom 
is not changed, play two important roles in the atomic cascade. First, the 
Stark mixing is nothing but the elastic scattering as long as only the rates of 
transitions between the /-sublevels of the same n are concerned. Second, the 
energy losses in elastic collisions lead to the deceleration of the epithermal 
exotic atoms. 

The basic features of the elastic scattering can be explained in the semi- 
classical approximation |23| . The motion of the [ip atom with the parabolic 
quantum numbers (ni,n2,m) in the electric field of a hydrogen atom is de- 
scribed by the effective potential^ 

where A = {rii — n 2 ), m w is the /ip reduced mass, C(r) = (1 + 2r + 2r 2 )e~ 2r 
is the electron screening factor, a is the electron Bohr radius. Neglecting the 
electron screening, this corresponds to the well known R~ 2 potential with a 
large dimensionless coupling constant 2Mg = ^nAM/m^ p ~ Mn 2 /m^ p ^> 1 
(M is the reduced mass of the fip + p system). The R~ 2 behaviour leads to 
the E~ x dependence of the elastic cross section. In the limit 2Mg ^> 1, many 
partial waves contribute to the scattering, and the differential cross section 
can be approximated by the classical result (g > 0): 

^ - f < 7 > 

, . n(n — 6) ^ 1 

w[ ~ ' = (2tt -ey e 2 sine ~ 32 sin 3 6/2 ^ 

The average coupling constant for given n is g = (n 2 — l)/2m Atp . The cross 
section (^^) has an unphysical singularity at 9 — > which must be regular- 
ized by taking the electron screening into account. 

3 The quantization axis is along the electric field, n = ri\ + 712 + \m\ + 1. 
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A semiclassical model of the Stark mixing was developed in and used 
with some variations and refinements in [|6|^,[T0|1. In this model, the exotic 
atom moves along a straight line trajectory with a constant velocity through 
the electric field of a hydrogen atom, and the Stark transitions induced by this 
electric field are calculated. Instead of straight line trajectories the classical 
trajectories can be used as in [PSR . The characteristic scale of the Stark 
mixing cross section is determined by the electron screening of the proton 
electric field, i.e. by the size of the hydrogen atom which makes Stark mixing 
the fastest collision process in the atomic cascade, see Figs. |[ 

Recently the problem of the elastic scattering of exotic atoms in excited 
states was treated in a quantum mechanical framework using the adiabatic 
expansion |^,^l|,|r5[, the calculation being done in a single-channel adiabatic 
approximation. The results show fair agreement with the semiclassical cal- 
culations [T^] for the collision energies E > 1 eV. A detailed discussion of 



the Stark collisions, including the relative role of transitions with different 
Al = If — li can be found in [ 11 ] . It would be desirable to extend the quan- 



tum mechanical model of the elastic collisions by including the effects of the 
shift and width of the nS states as well as the coupling between different 
adiabatic terms. 

The elastic scattering leads to the deceleration of exotic atoms if the ki- 
netic energy is larger than the target temperature. The effective deceleration 
rate for the \ip + p collisions is defined by the transport cross section: 

&n{E) = J (1 - cos 6)da n (E, 6) (10) 

where v is the atomic velocity, M MP and Mh are the [ip and hydrogen masses 
correspondingly. Using the classical result (|7],|8|) one gets the following esti- 
mationf] of the transport cross section: 

<(E) „ . (id 

The deceleration rate (Fig. 0) is a crude measure of the slowing down; in 
general one uses the differential cross sections in detailed kinetics calculations. 
The first estimation of the deceleration rates was done in ||12| , semiclassical 



l Note that the transport cross section is finite even without the electron screening. 
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calculations were done in [13, II |, and quantum mechanical calculations in 

BEL1IH- 



6 Coulomb De-excitation 



The Coulomb de-excitation (Table 1), where the transition energy is shared 
between the colliding particles, is an important acceleration mechanism pro- 
ducing 'hot' (E ^> 1 eV) exotic atoms (see Fig. |]a). The highly-energetic 
n~p atoms were discovered experimentally in p2]| , and a multicomponent 
structure of the energy distribution was found in [46]. Recent experiments 



in liquid and gaseous hydrogen |^7|,Q found the shape of the kinetic energy 
distribution to be consistent with the Coulomb mechanism. The cascade cal- 
culations show that the n~p atoms with kinetic energy E > 50 eV are 
not significantly decelerated between the Coulomb de-excitation and the nu- 
clear reaction. For the atoms with E < 20 eV the deceleration is important 
and the Coulomb peaks are expected to be smeared out. 

The earlier calculations of the Coulomb transitions []T6|-[T9|[ were rather 
controversial. The latest calculations [p0p|S2"| seem to clarify the theoretical 
picture. This removes one of the uncertainties in the cascade calculation of 
the past where a Coulomb rate-tuning parameter was used to normalize the 
calculated high-energy component to the measurements in pionic hydrogen 
|2D[ . It is not excluded that the Coulomb de-excitation can be part of a 
multistep process. For instance, the formation of excited molecular states 
can be followed by a Coulomb-like decay 0,|54|,^ . 



7 Atomic Cascade in H-D Mixtures 

The atomic cascade in mixtures of hydrogen isotopes provides additional 
information on the collisional processes in excited states. The muon transfer 
in excited states was studied in detail in the recent experiments [^D|,^T],^3| by 
measuring the relative yields of the up and \xd if-lines in liquid and gaseous 
HD mixtures. Figure [| shows a comparison of the experimental data with 
the cascade calculations for the probability Qis that the muon captured 
initially by the hydrogen in a HD mixture reaches the /ip ground state (the 
fraction (1 — Qis) is transferred to fid during the atomic cascade). Since 
the rates of the muon transfer (fip) n + d <-> (/irf) n + p are strongly energy 
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Figure 5: The Qis factor for the fip cascade in HD mixture vs. deuterium 
fraction Cd- The data are from |^0|,Q; the theoretical curves were 
calculated with the transfer rates |25| scaled by the factor k tr = 0.5 (solid 
lines) and k tr = 1 (dashed lines). 



dependent [g^-gg, the Qis factor is very sensitive to the kinetic energy 
distribution in the excited states p.p^.[26|, p7| . |3T l. 

The theoretical models tend to predict a stronger dependence of Qis on 
the deuterium fraction and the density than experimentally observed (this 
so-called "Qis-prcblem" is known since the first experimental results were 
obtained from the kinetic analysis of muon catalyzed fusion [[0], [EJ). The 
agreement between theory and experiment improves if the transfer rates are 
scaled by a factor of about 0.5 |BT[. This suggests that some mechanism may 



still be missing in the cascade model. One candidate is the resonant molecular 
formation in excited states (resonance side-path mechanism [[55|0 which can 
produce a significant inverse transfer /id — > fxp leading to an enhancement of 



Qis as discussed in [56 



8 Conclusion 

The current cascade models are advancing from a phenomenological level 
towards straightforward detailed calculations of the kinetics. The recent 
progress in quantum mechanical calculations of the scattering of exotic atoms 
in excited states makes it possible to get rid of tuning parameters and ad hoc 
assumptions in future calculations. The ultimate goal of ah initio cascade 
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calculations does not seem to be unrealistic any more. 

A detailed knowledge of atomic cascade is essential for precision X-ray 
spectroscopy of the 7i~p atom where the Doppler-broadening corrections must 
be taken into account in order to measure the nuclear width of the IS state. 



The new generation of precise experiments with the \ip and n~p atoms f49 
will make it possible to study the kinetics of the atomic cascade in great 
detail and to perform critical tests of the theoretical methods used in the 
calculations of various reactions with excited exotic atoms. The results of 
theses studies will also be very useful for other systems like K~p and pp. 
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